From canonical to nonautonomous solitons 



00 

o 
o 

(N 



oo 



Dun Zhao, 1 ' 2 Xu-Gang He, 1 and Hong-Gang Luo 2,3,4 

School of Mathematics and Statistics, Lanzhou University, Lanzhou 730000, China 

2 Center for Interdisciplinary Studies, Lanzhou University, Lanzhou 130000, China 

3 Key Laboratory for Magnetism and Magnetic materials of the 

Ministry of Education, Lanzhou University, Lanzhou 730000, China 

^Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100080, China 

In this paper we show a systematical method to obtain exact nonautonomous soliton solutions 
of the nonautonomous nonlinear Schrodinger (NLS) equation named by Serkin, Hasegawa, and 
Belyaeva [Phys. Rev. Lett. 98, 074102 (2007)]. We first perform the Painleve analysis of the 
nonautonomous NLS equation and obtain a compatibility condition. Under this condition a general 
transformation is then found, which can reduce the nonautonomous NLS equation to the standard 
one. From this transformation, all solutions of the standard NLS equation can be converted into 
the corresponding solutions of the nonautonomous NLS equation. This builds a unified picture of 
the nonautonomous solitons and the canonical ones and deepens the understanding of the nonau- 
tonomous soliton dynamics. Finally, the fundamental bright and dark solitons are taken as examples 
to demonstrate explicitly its applications. 
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Introduction. - The solitary wave is a fundamentally 
nonlinear object in nature. It has been widely observed 
in various fields since 1834 [1]. The concept of soliton 
was first proposed by Zabusky and Kruskal [2] since they 
found in 1965 that the nonlinear solitary waves can pass 
through each other and preserve their shapes and speed 
after a collision, just like an elementary particle. This 
special feature of these solitary waves explained the re- 
currence phenomenon observed by Fermi, Pasta, Ulam, 
and Tsingou [3] . Since then, extensive theoretical and ex- 
perimental investigations on the solitons have been made 
and some important applications of the solitons on the 
real world like the optical soliton communication [4] were 
also explored. 

Why can the soliton behave like an elementary parti- 
cle? As early as in 1870's, Boussinesq [5] and Raylcigh 
[6] have realized that the solitary water wave form results 
from the balance between the increase in local wave ve- 
locity associated with finite amplitude and the decrease 
associated with dispersion. In fact, this observation is 
also true to all solitons observed in various fields, namely, 
a localized nonlinear wave propagating in a medium can 
form a soliton when the effect of dispersion and that of 
nonlincarity is balanced. In addition, a dynamical bal- 
ance between dissipation and nonlincarity can also sta- 
bilize a solitary wave [7]. Therefore, the general picture 
that the balance between different effects in a nonlinear 
system can result in a solitary wave provides a good op- 
portunity to develop the concept of solitons in various 
complex situations. 

As the first extension of the soliton concept, in 1976 
Chen and Liu [8] found that the soliton can be acceler- 
ated in a linearly inhomogeneous plasma. At the same 
time the exact soliton solutions for the Korteweg-de Vries 
equation with varying nonlincarity and dispersion have 
also been found by Calogero and Degaspcris [9]. Re- 
cently, the soliton dynamics for the discrete [10, 11] and 



the generalized NLS equation with varying dispersion, 
nonlincarity, and dissipation or gain [12, 13, 14, 15, 16] 
have been extensively investigated. These solitons are 
apparently different from the classical solitons intro- 
duced by Zabusky and Kruskal [2] , which progagate with- 
out energy dissipation or gain [17]. In 2007, Serkin, 
Hasegawa, and Belyaeva [18] introduced a concept of 
nonautonomous solitons, which are the exact solutions 
of the nonautonomous NLS equation with time- and/or 
space-dependent dispersion and nonlincarity. Distin- 
guishing from the nonautonomous soliton, the classical 
one can be called as the canonical solition. The vari- 
ous complex nonlinear optical media [4] and the Bosc- 
Einstcin condensates (BEC) tuned by the Feshbach res- 
onances [19] provide many good grounds to study the 
dynamics of the nonautonomous solitons. 

Two questions arise: what is the essential property 
of the nonautonomous solitons? What is the relation- 
ship between the nonautonomous solitons and the canon- 
ical ones? Some preliminary studies based on Lax 
pair [18] and similarity transformation [20] techniques 
indicated that under some integrability conditions the 
nonautonomous solitons remain the basic property of the 
canonical solitons but both amplitudes and speeds of the 
solitons vary with time and space. However, the second 
question is not yet clearly answered. 

In this paper we focus on the second question and 
take a generalized NLS (or Gross-Pitaeskii) equation with 
varying dispersion, nonlincarity and harmonic external 
potential as an example. Using the Painleve analysis 
we first study a compatibility condition under which the 
nonautonomous NLS equation might be completely inte- 
grablc. Then, under this condition we try to construct 
a transformation which can reduce the nonautonomous 
NLS equation to a standard (or autonomous) NLS equa- 
tion. As a result, this transformation builds a mapping 
between the canonical solitons and the nonautonomous 



ones, which provides a systematical and significant way 
to study the dynamics of the nonautonomous solitons 
since the dynamics of the canonical solitons have been 
well-studied in more than hundred years. 

Model and Painleve analysis. - We start from a one- 
dimensional nonautonomous NLS equation 

du(x.t) ,, .d 2 u(x,t) . , ... . l9 . . 
i \: , +ef(x,t) fl %' ; + 6g(x,t)\u(x,t)\ 2 u(x,t) 
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+V(t)x 2 u(x,t) = 0, (1) 



where f(x,t) and g(x,t) denote the time- and space- 
dependent dispersion and nonlinearity, respectively, e 
and 5 are constants which describe different systems to 
which the nonautonomous NLS equation can be applied. 
V(t) represents a time-dependent confining harmonic ex- 
ternal potential strength. These coefficients are usually 
assumed to be real. 

Motivated by the connection between complete inte- 
grability and the Painleve property of partial differen- 
tial equations [21, 22], we perform the WTC test [21] to 
study possible integrability condition of Eq. (1). Follow- 
ing the Kruskal ansatz [23], u(x,t), v(x,t), f(x,t), and 
g(x, t) can be expanded on a non-characteristic singular- 
ity manifold tp(x,t) = x + <fi(t), where v(x, t) = u*(x,t) is 
viewed as an independent function. The expansions read 



u(x : t) = {x + 4>{t))-P Y. u m (t)(x + cj>(t)) m , 

"S° (2) 

V{x,t) = {x + ( j>{t))-1 E «m(*)(x + 0(*)r. 
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where uq e^ 0, vq ^ are required and p and q are integers 
to be determined. Likewise, 
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g(x,t) = E9n(t)(x + <j>(t)) n , 
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where /„(£) = ^ - ^' T> \ x =-4(t) and 9n(t) is similar. 

Inserting the above expansions into Eq. (1) and its 
complex conjugate equation, the leading order analysis 
leads to p = q = 1 and 2/o(i) + go(t)uo(t)vo(t) = 0. The 
recursion relations read 



30^0 



9ou"o 






(4) 



where Q m = (m — l)(m— 2)/o + 2go'«oUo an< ^ ^^ c expres- 
sions of F m (G m ) are shown in Ref. [24]. The zeros of the 
matrix determinant in Eq. (4) give the resonance points 
to = -1,0,3,4 of Eq. (1). While to = -1 and lead to 
trivial results, the compatibility conditions under which 
Eq. (4) has a nontrivial solution are given by 



m = 3 : vo (t) F 3 - u (t) G 3 = 0, 
m = 4 : wo (*) F 4 + u (*) G 4 = 0. 



(5) 
(6) 



These two conditions require that f(x,t) and g(x,t) are 
space-independent, i.e., f(x,t) = f(t) and g(x,t) = g{t). 
Furthermore, it is shown that Eq. (6) can be written as 



9 9 1 J 9 J 



0. 



(7) 



where the subscripts denote the time derivatives. This is 
a compatibility condition that Eq. (1) is Painleve inte- 
grable. It is interesting to note that this condition is com- 
pletely consistent with that obtained by Lax pair [18], 
but the space independence of dispersion and nonlinear- 
ity coefficients are required by the Painleve analysis. 

The complete integrability of Eq. (1) under the com- 
patibility condition can be further confirmed if there ex- 
ists a transformation that reduces Eq. (1) to the standard 
NLS equation 



i) 







—Q(X,T)+e—,Q(X,T) + 5\Q(X,T)\'Q(X,T) = 0. 



dT 



dX 



(8) 

Below we look for such a transformation in a general 
form of [25] 



u(x,t) = Q(X(x,i),T(i))e M(a; ' t)+c(t) , 



(9) 



where X(x, t), T(t), a(x, t) and c(t) are real functions and 
u(x,t) and Q(X,T) are required to be the solutions of 
Eqs. (1) and (8), respectively. Inserting Eq. (9) into 
Eq. (1) and comparing with Eq. (8), we obtain a set 
of differential equations, which have solutions under the 
condition Eq. (7), 



a{x,t) 
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4e/(t) \dt g(t)J VW 

-C\e J g j^dt' + C 2 , (10) 

x^-mf-^'fm*- (11) 

(12) 
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where C\ , C% , and C3 are arbitrary constants and are 
set to be zero in the following discussions. The Painleve 
integrability condition Eq. (7) together with the trans- 
formation Eqs. (9 - 13) consist of our central results of 
the present work. 

An important consequence of these equations is that it 
provides a systematical way to find the exact solutions of 
the nonautonomous NLS equation (1), as shown in Fig. 
1. For a given nonautonomous NLS equation, we first 
check if the coefficients satisfy the compatibility condi- 
tion Eq. (7). If it is true, then the nonautonomous NLS 
equation can be reduced to the standard NLS equation 
(8). All exact solutions, including the canonical solitons, 
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FIG. 1: Schematic diagram to obtain the nonautonomous soli- 
ton solutions of the nonautonomous NLS equation from the 
canonical solitons of the autonomous NLS equation. 



of the standard NLS equation can thus be converted into 
the corresponding solutions of the nonautonomous NLS 
equation. In this sense, a canonical soliton can be viewed 
as a "seed" of the corresponding nonautonomous solitons 
under the compatibility condition Eq. (7) , as emphasized 
on the title of the paper. 

Some remarks are in order, i) If f(t) = g(t) and 
V(t) = 0, the nonautonomous NLS equation Eq. (1) 
have the canonical soliton solutions (up to a phase) re- 
gardless of the explicit form of the time-dependent non- 
linearity and dispersion. This is because in this case the 
balance between nonlincarity and dispersion is not bro- 
ken down. In this sense the nonautonomous soliton is a 
quasi-canonical soliton. ii) When g(t) ^ f{t), the origi- 
nal balance between nonlincarity and dispersion is broken 
down. In this case the canonical soliton deforms itself to 
build new balance between nonlincarity and dispersion. 
In this sense, the nonautonomous soliton is a deformed 
canonical soliton. The amplitude of the soliton will be 
scaled by the factor of \/g(t)/f(t), as shown by c(t). iii) 
It is very interesting to note that the confining harmonic 
external potential is absent in the transformation equa- 
tions. However, the presence of the potential affects the 
balance between nonlinearity and dispersion and builds 
a deep connection between the optical solitons and the 
matter-wave ones. Once a BEC soliton is formed in a 
confining potential, it can behave like an optical soliton, 
possibly can be controlled by appropriate nonlincarity 
and dispersion managements even without the external 
potential, iv) If V(t) = 0, the solitons can be quasi- 
canonical or deformed depending on if fit) is equal to 
g{t) or not. On the contrary, if V(t) ^ 0, /(£) must 
be unequal to g{t). This leads to an important obser- 
vation that there does not exist the canonical and even 
the quasi-canonical the matter-wave solitons under the 
compatibility condition Eq. (7) and will be discussed in 
detail elsewhere. 

It is also helpful to mention some techniques to find the 
nonautonomous soliton solutions of the nonautonomous 
NLS equation in the literature. The Lax pair analy- 
sis is very useful in discussing integrability conditions 
[12, 16, 18]. The similarity transformation is usually lim- 
ited to the very explicit transformation form, as shown 



in Ref. [26]. Another similarity transformation reduced 
the nonautomonous NLS equation to a stationary NLS 
one has also been introduced [20]. Quite different from 
these techniques, the present work builds a connection 
between the nonautonomous NLS equation and its au- 
tonomous counterpart and provides a more systematical 
way to find solutions of the nonautonomous NLS equa- 
tion. To explicitly show how the method works, below 
we consider the fundamental bright and dark soliton so- 
lutions of the nonautonomous NLS equation without the 
harmonic external potential. 

Examples. - In this case, Eq. (7) has a solution 



f(t)=g{t)exp[-a g{t')dt'\, 



(14) 



where a is a constant and the transformation equations 
(10)-(13) become a(x,£) = -f exp(G Q (t))x 2 , X(x,t) = 
exp(G a (t))x,T(t) = fdt'g(t')exp(G a (t')), and c(t) = 
(l/2)G a (t), where G a (t) = af*g(t') dt'. 

When e = 1/2 and S = 1 (Se > 0), Eq. (8) 
has the fundamental canonical bright soliton solution 
Q(X,T) = sech(A") exp(iT/2) and when e = -1/2 and 
<5 = 1 (5e < 0), the fundamental dark soliton solution 
of Eq. (8) has the form of Q{X,T) = tanh(AT) cxp(iT). 
Starting from these two solutions, we show the corre- 
sponding nonautonomous soliton solutions of Eq. (1) 
for four different cases: g(t) — l,exp(i),exp(— t), and 
cos(i), which represent constant, enhancement, suppres- 
sion, and periodic modulations of nonlinearity, respec- 
tively. The corresponding dispersion modulations follow 
Eq. (14). It is noted that a = leads to G (i) = 0, 
which is trivial up to a phase, as mentioned above. A 
nonzero a has nontrivial results and without loss of gen- 
erality we take a = 1 below. In Fig. 2 and Fig. 3 we 
explicitly present the nonautonomous bright solitons and 
the nonautonomous dark solitons for four different non- 
lincarity modulations, respectively For comparison, we 
also plot the canonical solitons, as shown in Fig. 2(a) 
and Fig. 3(a), respectively. 

Fig. 2(b) and (c) show that the canonical bright soli- 
ton becomes more and more either sharper or broader 
depending on enhancement or suppression of the non- 
linearity. When the nonlinearity keeps unchange and the 
dispersion is suppressed, the canonical bright solition also 
becomes more and more sharper, as shown in Fig. 2(e). 
This can be understood by the fact that the soliton is 
due to the balance between dispersion and nonlinearity. 
Most interesting case is Fig. 2(d), where the nonlincar- 
ity modulation is periodic. As a result, the canonical 
bright soliton is also modulated periodically These all 
results indicate that the nonautonomous bright soliton 
and its canonical counterpart has a close relationship. 
In all cases, it is easy to check that the integrations of 
J \u(x,t)\ 2 dx keep unchange in time. The similar anal- 
ysis is also true to the nonautonomous dark solitons, as 
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FIG. 2: (a) The canonical bright soliton of Eq. (8) with e = 
1/2 and 8—1 and the corresponding nonautonomous bright 
solitons of Eq. (1) for different nonlinearity modulations: (b) 
g(t) = exp(i), (c) g(t) = exp(-i), (d) g(t) = cos(f), and (e) 
g(t) = 1. In all cases a = 1. 



However, the presence of dissipation and/or gain in these 
systems will complicate the present analysis and a non- 
trivial extension is needed. 

Finally, it should be pointed out that the present anal- 
ysis can also be applied to all exact solutions of Eq. (8), 
including the multi-soliton cases. This provide a system- 
atical way to study the dynamics of the nonautonomous 
NLS equation. 

In conclusion, based on the Painlcvc analysis and a 
general transformation obtained we have reduced the 
nonautonomous NLS equation to the standard one. This 
provide a systematical and significant way to obtain the 
nonautonomous soliton solutions of the nonautonomous 
NLS equation from its autonomous counterpart, i.e., 
the canonical solitons. The nonautonomous solitons ob- 
tained result from a balance between dispersion, nonlin- 
earity, and/or an external potential applied, just like the 
canonical solitons. This result builds a unified picture of 
the nonautonomous and the canonical solitons and pro- 
vides an effective way to control the nonautonomous soli- 
ton dynamics. 
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FIG. 3: The canonical dark soliton of Eq. (8) with e = -1/2 
and S — 1 and the corresponding nonautonomous dark soli- 
tons of Eq. (1) for different nonlinearity modulations same 
as those in Fig. 2. For clarity, the amplitude of the nonau- 
tonomous dark solitons is normalized by r\t = exp \\G\[tM. 



shown in Fig. 3. These results shed light on the un- 
derstanding of the nonautonomous soliton dynamics and 
provide an exact way to make a dispersion and/or non- 
linearity management of solitons. It is expected to have 
a realistic application to the optical soliton communica- 
tion technologies and the matter- wave soliton dynamics. 
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